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We predict that the mesoscopic tensile force fluctuations 
in metal quantum point contacts (nanowires) arise as a result 
of finite electric voltage on the contact. They are due to re- 
configuration of the electronic subsystem and are correlated 
with the nonlinearities of the current-voltage characteristics 
of the contact. The observation of the effect would directly 
confirm the recently suggested "free-electron" mechanism of 
mesoscopic force fluctuations observed in nanowires under de- 
formation. 

The related magnetic susceptibility fluctuations and role of 
topology of the wire cross section are discussed as well. 

The quantization of electrical conductance of 3D quan- 
tum point contacts (QPC) jjj has been observed in a 
variety of metal mesoscopic contacts @. The simplest 
model of the 3D QPC is a "nanowire" of length L and 
diameter d ~ \p, connecting two bulk conducting reser- 
voirs, where A^ is the Fermi wavelength in the system 
(see Fig. |l|), and the underlying atomic structure is the 
positively charge "jcllium". It is quantized character of 
transverse motion of electrons that is revealed in conduc- 
tance quantization as a function of d. Since the applied 
driving voltage changes the population of occupied sub- 
bands, nonlinear current - voltage dependence quantum 
point contacts, was predicted j^-ll. The theory is in a 
good agreement with experiment on 2D and 3D systems 
0. On the other hand, the nonlinear conductivity of 3D 
bismuth QPC || shows only qualitative agreement with 
the theory, while measurements on 3D gold QPCs || did 
not show the predicted type of nonlinearities altogether 
©■ 

The experiments on metal QPCs under deformation 
showed that the mechanical stress in the wire fluctuates 
as a function of its elongation, and the fluctuations are 
correlated with the conductance jumps Q. This seemed 
to require a more complex mechanism, and the often- 
used explanation of these phenomena invokes the atomic 
rearrangement processes, and is supported by molecular 
dynamics calculations (e.g., fll^l ). 

Recently, an elegant alternative explanation of the 
force fluctuations was suggested, based on the reaction 
of the free electrons to the mechanical deformation of 
the contact region |l3| . The role of the atomic structure 
of the wire was again reduced to providing a "jellium" 
background. In this "free-electron" model, the longitu- 
dinal force, being the coordinate derivative of the ther- 
modynamic potential, is sensitive to the positions and 
occupancy of electronic subbands in the wire. The lat- 
ter depends on the shape of the cross section and on the 



elongation of the wire (assumed to take place at a con- 
stant volume). The positions of cusps in F as a function 
of the elongation would naturally coincide with those of 
the conductance steps. 

In this paper we show that the "free-electron" mecha- 
nism of tensile force fluctuations will lead to related ef- 
fects: mechanical force and magnetic susceptibility fluc- 
tuations in a nanowire as a function of applied driv- 
ing voltage at the same values of eU, as the features 
of the differential conductance, Gdis{eU). Investigation 
of these effects can be done using existing experimental 
techniques and would provide an independent test of the 
mechanism suggested in [ jl3"| , and confirm the decisive 
role of electronic subsystem in determining both trans- 
port and mechanic properties of metal quantum contacts. 

We start from the expression for the grand potential of 
the electronic subsystem at zero temperature and voltage 
for a wire of uniform cross section (Blom et al. 

n 

x6(E F - E n (L,V)). (1) 

Here m* is the effective mass of an electron, 6(x) is the 
Heaviside step function, and E n is the energy of nth elec- 
tronic transverse mode in the wire, which depends on 
the wire length, L, and volume, V. We assume that the 
length of the wire is much larger than its diameter d, 
which allows us to set the electrical potential to zero in 
the wire. (Due to screening in metal the effects of finite 
bias will be felt only at distances d from the ends, see 
H and references therein.) 

For the sake of simplicity, we will neglect both elastic 
and inelastic scattering in the nanowire. Their contribu- 
tions are of order (L,d)/l s , where l s is the corresdponding 
scattering length The condition (L, d)/l s <C 1 is sat- 
isfied both for the elastic and electron-phonon scattering. 
In the latter case we use the estimate l e -ph ~ hv F / Xlod, 
where u)p> is Debye frequency, and A < 1 is the electron- 
phonon coupling constant. The electron-electron scatter- 
ing length can be estimated as l e -e ~ TiVp /ep (ep /eV) 2 , 
the bias eV playing the role of effective temperature. 
For the effects of electron-electron scattering to be small 
we need d/l e - e <C 1 (because the longitudinal momen- 
tum conservation and transverse quantization in the elec- 
tronic subsystem suppress the electron-electron scatter- 
ing inside the nanowire). The corresponding restric- 
tion on the applied voltage is eV < ep/s/Nx, where 
N± ~ d/Xp is the number of transverse channels in the 
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quantum contact. This condition is compatible with bias 
being of the same order as the interlevel spacing in the 
contact, Ae ~ e F /N±, which is necessary for the obser- 
vation of nonlinear effects discussed in this paper. 

Under the above assumptions, we can consider right- 
moving and left-moving electrons as two independent 
subsystems, with chemical potentials /il and fiR respec- 
tively (Fig. @): 
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(1 - P)eU; /j, r = E f + (3eU. 
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The grand potential thus becomes 
n eU = -(Cl(E F + f3ell) + n(E F 



(l-p)eU)). (3) 



Parameter determines the asymmetry of the voltage 
drop on the contact. (Usually symmetric voltage drop is 
assumed (J3 — 1/2), in which case the differential con- 
ductance is always a multiple of \Gq. || Generally 
can deviate from 1/2 ( 0) and be voltage-dependent.) 
It should be in principle determined self-consistently by 
solving the electrostatical problem for the wire and its 
surroundings at given eU JhJ. We consider here two lim- 
iting cases: (a) symmetric voltage drop, = 1/2, and (b) 
perfect screening. In the latter case /3(eU) is determined 
from the condition of no charge accumulation in the wire, 

N(E F + 0eU) + N(Ep - (1 - 0)eU) = 2JV(0), (4) 

where 



n(e) = 2lj ^Y. e ( E - E ^y E ~ ( 5 ) 



It is easy to see, that the differential conductance of 
the system is given by 

AT A 2p r EF+peU ^ 
G diS (eU) = £- = -±f dE-£e(E-E n ) 

dU dU h J EF -(i-p)eU „ 



Gq ^{0O(Ef + /3eU - E n ) 



Here E — n 2 h /(2m*d ), and j mn is denotes the nth 
positive zero of the Bessel function J m (z); m is the mag- 
netic quantum number. The mechanical force is given 
by 



F(u) = - 
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dL 
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+f(e F -(l-0(u))u;e mn )}, (8) 



where e F = E F /E ; u = eU/E ; e mn = 
E mn /E ; F a = 7r 2 h 2 /(2m*d 3 ) 1 and f{x;y) = 
((2/3)(x - y) 3/2 - (x - y) 1/2 y) 6(x-y). The nonlinear 
dependence F(eU) and Gdis(eU) in both limiting cases 
is shown in Fig. |^. The nonlinear conductance is strongly 
dependent on the character of screening in the wire. On 
the other hand, the qualitative character of the force 
fluctuations is the same, and mechanical force and dif- 
ferential conductance still show singularities at the same 
applied voltages in both limits Jig ]. The absolute mag- 
nitude of force fluctuations for d ~ lnm is of order 1 nN, 
in agreement with earlier results |ll],|l3| . 

Another way of changing positions of quantized lev- 
els, and thus the properties of the contact, is by ap- 
plying longitudinal magnetic field. The characteristic 
field sweep scale, corresponding to interlevel spacing, is 
though ~ Q /d 2 , where $o = hc/e is magnetic flux quan- 
tum. (l],|(|[rj| For a metal QPC with d ~ 1 nm this yields 
unrealistic fields of order 10 3 T. This means that in metal 
contacts, appreciable dependence of contact's properties 
on the magnetic field can take place only when the Fermi 
level is already very close to one of the quantized en- 
ergy levels. This can be achieved, e.g., by mechanical 
deformation of the contact, or by applying finite driv- 
ing voltage. We will concentrate on the latter possibility, 
which is reversible and promises better opportunities for 
the necessary fine tuning. 

The magnetic field can be taken into account in the 
perturbation theory Jl7| ] , valid in the limit of weak field 
(large cyclotron radius, r c 3> d) |l8|]one finds for the 
quantized transverse levels 
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+U -£{9(E F + 0eU - E n ) ~ 0(E F - (1 - 0)eU - E n ))}, 
du 

where Gq = 2e 2 /h is the unit quantum conductance. 
Therefore Gdis(eU) shows a structure at the volt- 
ages when consequent transverse energy levels enter the 
current-carrying interval [E F — (1 — 0)eU, E F + 0eU]. In 
the limit of perfect screening 0(eU) was found numeri- 
cally for both models we considered: the wire of square 
cross-section d x d, and the roun wire of diameter d. 
The quantized levels in the wire are given by 



Here r\ is the magnetic field measured in units of H = 
(hc/e)/{-Kd 2 /4). Since r\ <ti 1, we keep in @ only linear 
in i] terms, including the spin splitting (the last term in 
the parentheses). Here s = ±1 is the projection of spin, 
g is the g-factor, and mo is the free electron mass. 

The differential conductance and force fluctuations are 
thus given by 
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F(r,,u) = -F E E {f(tF + l3u-i s m M) (11) 

s— ±1 n— 1 m— — oo 

+/(e F -(l-^;e^ n (»j))}, 

where e mn = E mn /Eo. The factors of one half before 
Gq , Fq reflect the spin splitting in the magnetic field of 
previously degenerate energy levels. 
The magnetization of the wire is 
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MM = -v[dH VT = ^ B 22 

x {(e F + u/2 - i s mn ( V )) 1/2 8 (e F + u/2 - e^fa)) 
+ (e F - u/2 - i s mn (v)) 1/2 6 (e F - u/2 ~ e s mn ( V ))} , 

where fig = eh/(2m c) is the Bohr magneton. 

The effects of the applied weak magnetic field are 
described by magnetoconductance coefficient, cr(u) = 
{dG<mff/dH) VT . H _ , magnetotension coefficient, T(u) = 
(9F/dH)y T . H _ , and the magnetic susceptibility, 
x (u) = {dM/dH) VT . H=Q . 

Keeping only the singular terms, we find the following 
expressions: 



a(eU) = -\g q E 
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It demonstrates a series of inverse square root singular- 
ities at the same values of driving voltage, as the fea- 
tures of differential conductance and force fluctuations 
(see Fig. ||). The features of x( e ^0 are better pronounced 
than those of the former coefficients, which could out- 
weight the small magnitude of the effect and make mea- 
surements of magnetic susceptibility of a metal point con- 
a more sensitive tool for investigation of electronic 
density and potential redistribution in metal point con- 
tacts. 

In conclusion, using a simple model, we showed that 
finite driving voltage can lead to mechanical force fluc- 
tuations and singularities of magnetic susceptibility in 
metal quantum contacts. The mechanism of these ef- 
fects is voltage-induced nonequilibrium redistribution of 
electrons over quasi-lD subbands in the contact. On 
the other hand, magnetoconductance and magnetoten- 
sion coefficients are shown to be exactly zero, and the 
corresponding effects to be at least of order (H/Hq) 2 , 
where Hq m 10 3 T in a nanometer size contact. 
(13) Experimental investigation of the predicted effects 
would clarify the role played by electronic subsystem in 
behavior of metal quantum contacts. 

I am grateful to I. Affleck, E. Bogachek, A. Bratkovsky, 
and S. Rashkeev for helpful discussions. 
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Note that the magnetconductance and magnetotension 
coefficients contain the first power of de s mn {-q) /drj. There- 
fore they are exactly zero, due to cancellation of terms 
with opposite m, s: 

a = 0; T = 0. (15) 

This means, that the magnetoconductance and magne- 
totension in a metal quantum contact are the effects of 
second order in rj <C 1. They would thus appear as nu- 
merically small, extremely narrow peaks in voltage de- 
pendence of the corresponding functions, and cannot be 
in a satisfactory way investigated in our simple model, 
neglecting the effects of finite temperature and scatter- 
ing. 

On the contrary, the magnetic susceptibility contains 
the second power of de s mn {rj) / 'drj , and is thus nonzero: 



(14) 
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FIG. 1. Schematic view of a metal quantum point contact 
(nanowire) and distribution function of right- and left-moving 
electrons in the wire at finite driving voltage. Solid lines cor- 
respond to quantized values of transverse momentum. 




2 eU/E 4 



4 



FIG. 2. Differential conductance and force fluctuations 
vs. applied voltage in a nanowire of square cross section 
d x d (a,c) or circular cross section of diameter d (b,d). The 
force and bias are measured in units of Fq = -k 2 Ti 2 / '(2m* d 3 ) 
and Eo = n 2 h 2 /(2m*d 2 ) respectively. The Fermi energy is 
Ef — 11-Et)- Solid line: symmetric voltage drop (/3 = 0.5). 
Dots: perfect screening. 
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FIG. 3. Magnetic susceptibility vs. applied voltage in a 
nanowire of round cross section. The unit xo = m^B /(jn* Ho), 
where Ho = 8hc/(d 2 e). We chose g — 2, m* = m e . 
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